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Abstract 

Quantum Chromodynamics on a lattice with Wilson fermions and a chirally twisted 

Qh! mass term for two degenerate quark flavours is considered in the framework of chiral 

CD . 

r~ | . perturbation theory. The pion mass and decay constant are calculated in next-to-leading 

order including terms linear in the lattice spacing a. 

$_i ■ Numerical simulations of Quantum Chromodynamics with dynamical quarks are facing a dif- 



ficulty with light quark masses. Due to a severe slowing down of the common algorithms for 
small quark masses it is presently not possible to implement realistic values for the masses of 
the u- and d-quarks. In this context chiral perturbation theory [U El IS] is a very useful tool. It 
amounts to an expansion around the chiral limit where the quarks are massless. Meson masses 
and other physical quantities are then expanded in powers of quark masses (modified by log- 
arithms). In its range of applicability, chiral perturbation theory can be used to extrapolate 
the results of Monte Carlo simulations of QCD into the region of small quark masses. On the 
other hand, numerical simulations of lattice QCD can provide information about low energy 
parameters of chiral perturbation theory, the Gasser-Leutwyler coefficients, see |4| for a review. 

When chiral perturbation theory is applied to lattice QCD, the lattice spacing a enters as 
a second expansion parameter. Correspondingly, the low energy effective Lagrangean contains 
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additional terms proportional to powers of a, representing the lattice artifacts Physical 
quantities appear in double expansions in quark masses and in a. A comparison of numerical 
results with chiral perturbation theory in next-to-leading order has been made in 

An attractive framework for numerical simulations of QCD with Wilson fermions is the 
lattice formulation with a chirally twisted quark mass matrix [HI El- It allows to avoid problems 
with quark zero modes and simplifies the renormalization procedure. In view of Monte Carlo 
calculations in twisted mass lattice QCD it is desirable to extend chiral perturbation theory to 
this case. 

In this letter we consider lattice QCD with Nf = 2 quark flavours and a twisted mass term. 
For simplicity we restrict ourselves to the case of degenerate quark masses (m u = ma = m). 
The pion mass and the pion decay constant are calculated in chiral perturbation theory in 
next-to-leading (one-loop) order, including lattice terms linear in the lattice spacing a. A new 
feature of the twisted mass case is the fact that the minimum of the effective Lagrangean is 
shifted by an amount proportional to a and the expansion has to be rearranged correspondingly. 

In lattice QCD a chirally twisted mass term can be introduced in the form qM(uj)q with 

M(uj) —m'e 1W75T3 = rnf cos(u;) 1 + im' sin(u;) 75 r 3 = m + i ^75X3 , (1) 

where 

m — m! cos(o;) , \i = m'sin((j) , (2) 

and Tfc are Pauli matrices. In the continuum the chiral rotation of the mass term can be removed 
by a chiral transformation of the quark fields according to 

g = e -i-75T 3 /2 g / ; (3) 

which leaves the kinetic term invariant. On the lattice, however, it is not possible to remove 
the chiral twist due to the presence of lattice artifacts. 

Chiral perturbation theory is based on the low energy effective Lagrangean, which describes 
the dynamics of pions in terms of SU(2)-valued matrices 

U(x) = exp ir b (x)nj , (4) 

where 7r&(x) are the pion fields. In the continuum it is given to leading order by 

£2 = C 2Mn + C m = ^Tv (fytft &>U) - ^ Tr ( X C/ f + U X ] ) ■ (5) 
The chiral symmetry breaking mass term contains the matrix 

x = 2B M = 2B ml = X0 l, (6) 

with 

B = J_ <o|im|0) = i (o|dd|o) . (7) 



To lowest order in m the pion mass is given by 

- 2 - - - 2B m. 



Xo 



(8) 



Lattice QCD can be treated in chiral perturbation theory by adding terms which represent 
the lattice artifacts. The chiral Lagrangean to lowest order is extended to 

C 2 = C 2Mn + C rn -^Tv(pW + U^) , (9) 

where 

p = 2W a 1 (10) 

and the coefficient Wo parametrizes the lattice artifacts. 

Now let us consider the effective Lagrangean for QCD with a twisted mass term. The 
twisting of the mass term 

x _> x (u) = 2B e-' luJT3/2 Me~ iuJT3/2 = 2B me- iujT3 (11) 

can be undone in the kinetic and mass terms by a chiral transformation of the axial type 

U = e - iuJT3/2 U'e' iujTa/2 . (12) 
The lattice term, however, is transformed to 

Tt(p(u)U* + U'p{uy) (13) 

with the matrix 

p(u) = e^p. (14) 

In this way the twist can be transferred to the lattice term. Omitting the primes on the chiral 
field U, the leading order effective Lagrangean reads 



C 2 = S Tr farf &>U) - ^ Tr (( X + p(u)) E7+ + U( X + p(w))t) 



with 



X = 2B m'l. 

In next-to-leading order the following terms contribute 0: 

U = ^fTv(d fl Ud^)-^Tv(x^ + Ux^-^fTr(p(uj)U^ + Up(ujy 
-L x [Tr (dpUdpUt)] 2 - L 2 Tr (d^Ud^) Tr (dJJd v tf 
-U Tr 



(15) 
(16) 

(17) 



dpUdprf] J + U Tr {dJJdjfl) Tr ( X C/ f + U X ] ) 
+W 4 Tr (dpUdprf) Tr (p{uu)U ] + Up{wf) + L 5 Tr (d^Ud^ [ X U^ + U X ] 



+W 5 Tr (dJJdjfl [p(u)W + Up(uj] 
-W e Tr (xf/ f + U X ] ) Tr (p(cu)t/ t + Up{ 



— L 

f 



Tr ( X tft + f/x 1 



L 7 



Tr - f/ X ! 



-W 7 Tr (xU ] - U X ] ) Tr (p(c;)f/ t - Up( 
-W 8 Tr (x[/V(^)^ f + U P (uj) ] Ux ] ) + C(a 
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-cV) - 
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LgTr^x^ + ^f/x" 



Spurion analysis shows that the twisting of the mass term does not produce further terms in 
£4, but amounts to replacing p by piui). 

In contrast to the untwisted case, the minimum of the effective action is not located at 
vanishing fields 7r fo = 0, corresponding to U = 1, but at the point 



7T3 = ^0 ^ — : ssnu 
B m' 



l--^-(4L 6 + 2L 8 -2W 6 -W 8 ) 
^0 



+ 0(a 2 ), 7fi = 0, n 2 = 0. (18) 



Therefore a shift 

7T 3 = < + 7f 3 (19) 

has to be performed before expanding the effective action. 

Using the resulting expression we calculated the pion propagator in next-to-leading order. 
From it the pion mass and wave function renormalization are obtained. For the pion mass we 
get 

ml = Xo + Po + 8^m + 2L r 8 -2Ll-L r B ) 

+8^(4iy 6 r + 2W r 8 - 2W r A - W; - 2L\ - LI) 
(Xo + Pol j ( Xo + Po \ m 



where 



Xo = 2B m', p = 2W acosco , (21) 



L r k are the renormalized chiral parameters and A is the renormalization scale. The wavefunction 
renormalization in dimensional regularization is given by 



1 - ^2X0(2^4 + L B ) - (2W 4 + W s ) 



24vr 2 F r 2 



(Xq + Po) 



2 , 1 (a \ , 1 i fXo + Po 
- + ln(4-7r) - 7 + 1 - In I — — — 



(22) 



Setting the twist angle u> to zero, the expression for m 2 is consistent with the result of [Zj- 
Apart from numerical factors which depend on Nj, to this order of chiral perturbation theory 
the chiral twist amounts to adding the factor cosu to p. It should be noted, however, that 
this rule will not apply to higher orders, where the above mentioned shift in the pion fields 
introduces new vertices. 

In the case of maximal twist, u = ir/2, the lattice artifacts vanish to the order of the 
expansion considered here, as has been observed for lattice QCD in general in [Tfl| ITT]. 

The second physical quantity we calculated is the pion decay constant F n , given by 

(0|J£°|7r 6 (p)) =iF vI ?5 ab , (23) 
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where J a is the axial current. The Noether procedure yields an expression for the axial current, 
which depends on the coefficients and Wk and contains products of the pion fields. A one- 
loop calculation gives 

F* = Fo (l + ^[xo (2LI + LI) + p (2WJ + W£)} - 

In the limit a — > this coincides with the result of [B]. 

The formulae given here can be used in the analysis of numerical results from unquenched 
simulations of twisted mass lattice QCD with two light flavours. A useful technique in numer- 
ical simulations of QCD is partial quenching of quarks. In view of applications to partially 
quenched simulations it suggests itself to generalize the results above to partially quenched 
chiral perturbation theory [T2*] . 

We thank I. Montvay, E. Scholz and F. Farchioni for discussions. 



(Xo + Po) In 



Xo + Po 
A 2 



(24) 
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